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Abstract 

Recent lattice calculations performed at zero temperature and in the max- 
imal center gauge indicate that quark confinement can be understood in this 
gauge as due to fluctuations in the number of magnetic vortices piercing a 
given Wilson loop. This development has led to a revival of the vortex con- 
densation theory of confinement. For a SU(2) gauge group, we show that 
also at finite temperatures, center vortices are the relevant collective infrared 
degrees of freedom determining the long-range static quark potential; in par- 
ticular, their dynamics reflect the transition to the deconfining phase. 
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Introduction. One of the most intriguing prospects of strong interaction physics 
is the expected existence of a deconfined phase above a certain transition temper- 
ature. While strongly interacting matter, i.e. matter carrying a color quantum 
number, has to date only been observed in the form of color-singlet bound states 
called hadrons, the deconfined phase is characterized by the possibility of colored 
constituents propagating freely over distances large compared with hadronic sizes. 
Upcoming experiments at RHIC and LHC are expected to probe the deconfined 
regime, and thus efforts to establish an adequately detailed theoretical picture of 
the deconfinement phase transition and the deconfined phase acquire a measure of 
urgency. 

The investigation of the deconfinement transition is irrevocably tied to the question 
of the mechanism leading to confinement in the low-temperature phase of QCD, 
the gauge theory of the strong interactions. It should be noted that in different 
gauges, different characteristic infrared collective degrees of freedom dominate the 
physics of confinement, and the convenience of a gauge hinges on the ease with 
which these infrared degrees of freedom can be identified and isolated. Examples 
of gauges which are in this sense convenient are the maximal Abelian gauge, which 
displays monopole dominance and thus underpins the dual superconductor picture 
of confinement and the maximal center gauges [^-0, which display center 

dominance and lead to the center vortex picture of confinement It is this 

latter picture we concentrate on in this letter. 

Maximal center gauges are defined in the lattice formulation of Yang-Mills theory by 
the requirement to choose link variables on the lattice as close to center elements of 
the gauge group as the gauge freedom will allow. There are many slightly differing 
ways of implementing this general idea; in this work, the specific prescription called 
"direct maximal center gauge" was adopted, for details see e.g. 0. Furthermore, 
in what follows, we will be concerned with an SU(2) gauge group, the center of 
which consists of the elements ±1. The maximal center gauges have the following 
convenient property: If one replaces the gauge-fixed link variables with the center el- 
ements closest to them (so-called center projection), then a Monte Carlo calculation 
of the long-range zero-temperature static quark potential via the Wilson loop yields 
virtually the same result as when using the full unprojected link variables 0. This 
empirical observation is called center dominance. Note that in both cases, the weight 
in the Monte Carlo calculation is the full Yang-Mills action; only the observable, 
i.e. the Wilson loop, is sampled using either the full or the projected configura- 
tions. Center dominance is interpreted as meaning that the maximal center gauges 
successfully concentrate the information relevant for confinement on the center part 
of the link variables, thereby rendering the projection procedure of discarding the 
residual deviations from the center elements inconsequential. 



2 



In order to further physically interpret the center projected link configurations, one 
can connect them with a notion of vortices as follows. After center projection, 
one is left with lattice links associated with center elements ±1 of the gauge group 
SU(2). If the links bordering any given plaquette on the lattice multiply to — 1, 
then a vortex is said to pierce that plaquette. The vortices defined in this way form 
closed two-dimensional surfaces in a four-dimensional space-time lattice, or closed 
(magnetic fiux) lines on the three-dimensional lattice describing, say, one time slice. 
An important property of these vortices is that they contribute a factor —1 to any 
Wilson loop whenever they pierce its minimal area; in this way, they produce the 
same value for the Wilson loop as one obtains by multiplying the original center- 
projected links making up the loop (this is simply Stokes' theorem). If a Wilson 
loop area is pierced by vortices an even number of times, then it takes the value 
if a Wilson loop area is pierced by vortices an odd number of times, then it 
takes the value —1. Thus, if the distribution of vortices in space-time fiuctuates 
sufficiently randomly, then a strong cancellation will occur between configurations 
where Wilson loops are linked by an even or by an odd number of vortices; an ever 
better cancellation as the area of the Wilson loop increases leads to an area law decay 
of the Wilson loop expectation value, which implies a linear static quark potential. 
This is the essence of the random vortex picture of confinement. After having 



been proposed as early as 1978 [|1T|-[]13|, glimpses of such vortex configurations were 



afforded by the Copenhagen "spaghetti" vacuum and several investigations 
were devoted to identifying vortices on the lattice, leading to different definitions 
of vortices [|15|,|l^,0. The recent efforts in the framework of the maximal center 



gauges outlined above f^-^ have given new impetus to this particular description of 
confinement, foremost due to the verification of center dominance ||^. Subsequently, 
further properties of center vortices have been investigated, e.g. the connection 
to the monopoles of the maximal Abelian gauge 0, and the observation that the 
density and binary correlations of these vortices display the correct renormalization 
group scaling properties 0,0. This means that they survive the continuum limit 
and represent physical objects rather than lattice artifacts. 

Having identified such a viable picture of confinement, it is natural to confront it 
with the question of the deconfinement transition. In this letter, we verify that the 
center vortex picture makes sense also at finite temperature, since center dominance 
persists. Also the transition to the deconfined phase with a vanishing string tension 
is correctly reproduced. We further explore properties of vortices as the temperature 
is raised across the deconfinement transition and we cast the deconfinement order 
parameter into the vortex language. 

Zero temperature vortex properties. For later reference and in order to enable 
the reader to estimate the accuracy which is achieved in our numerical calcula- 
tions, we briefiy summarize our zero temperature results. Our numerical data were 
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Figure 1: Scaling of the planar density of vortex intersection points with a given 
space-time plane. 



obtained using lattices of size 12^ and 16^. The string tension is extracted from 
center-projected links by calculating nx n Creutz ratios, n = 3,4,5. One observes 
that the results for the Creutz ratios for different n agree already for these small n, 
indicating that the Coulomb part of the heavy quark potential is absent when the 
latter is evaluated using center-projected links |Q. The Creutz ratios allow us to 
extract (ja?{j3 = 2.3) = 0.12(3). The /3-dependence of aa^ is compatible with one- 
loop perturbative scaling. Two-loop effects are small compared with the statistical 
error of our numerical data. Furthermore, for the lattices used (12^ and 16^), the 
finite size dependence of dimensionless quantities was negligible compared with the 
statistical errors. However, as is well known, considerable finite size effects persist in 
the extraction of the mass scale, as e.g. encoded in the lattice spacing, as a function 
of (3. In the large scale analysis [|T^], an interpolation of the numerical data with the 
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help of one-loop scaling to the renormalization point (3 = 2.3 reveals aa"^ = 0.136 
for a 10^ lattice, cxa^ = 0.121 for a 16^ lattice, while aa^ = 0.107 is found for a 32^ 
lattice. Compatible with this, we use throughout this letter aa?{[3 = 2.3) = 0.12 and 
a = (440 MeV)^ as reference scale for assigning physical units to measured quanti- 
ties. When finite size uncertainties are mentioned below, this always refers to the 
aforementioned difficulty in defining the mass scale at a given (3. 

In order to extract vortex properties, we fix the configurations to the direct maximal 
center gauge by using the procedure described in |]^. To ehminate Gribov copies, 
we perform random gauge transformations and repeat minimizing the gauge fixing 
functional. It turns out that three of these runs are sufficient at zero temperature 
For temperatures close to the transition point, we observe an increasing influence of 
the Gribov copies making six gauge fixing iterations necessary. The scale dependence 
of the vortex density is shown in figure |l|. We find a perfect scaling of the planar 
density p of vortex intersection points with a given space-time plane in the range 
/3 e [2.25, 2.55] for lattice sizes 12^ 16^ We finally extract p = 3.6 ± 0.2fm-l The 
large error bar in p results from the uncertainty in the magnitude of finite size effects. 
The above value, which is larger than the one quoted in |^ , is now in agreement with 
the one quoted in within the error bars. The faulty estimate for p given in 
is due to the fact that the Coulomb part of the full zero temperature static quark 
potential was underestimated in which subsequently led to an overestimate of 
the reference scale. As a consequence, mass scales in and also in should 
be rescaled upwards by a factor 1.35. Note also that the new corrected value of 
cr/p = 1.4 ± 0.1 as opposed to the old value cr/p = 2.5 implies that a model of 
random intersection points of vortices with a two-dimensional space-time plane 
overestimates, rather than underestimates, the string tension as a function of the 
planar density of points p. Such a model leads to a value of a/p = 2. However, 
this does not affect the motivation for the lattice measurements in ^ nor their 
validity, up to the aforementioned trivial rescaling of mass scales. In ||^, correlations 
of an attractive type were observed between intersection points of vortices with a 
two-dimensional space-time plane. It seems plausible that such correlations curtail 
the randomness of the distribution of the points and thus reduce the string tension. 
Indeed, further below it will become clear that it is a pairing of the vortex intersection 
points which ultimately leads to deconfinement. However, one should bear in mind 
that the connection between correlators such as measured in |^ and the string 
tension is very indirect and one can not in general predict the behavior of the string 
tension from these correlators alone. 

Center dominance at finite temperatures. Finite temperature Yang-Mills the- 
ory can be cast into a path integral formulation on a space-time manifold of finite 
(Euclidean) time extension 1/T, where T is the temperature. Gauge fields on the 
lattice obey periodic boundary conditions in the time direction. The static quark 
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potential V{r) is calculated from the correlator of two Polyakov loops [0, 



{L{x)L{y)) oc exp{-V(r)/T}, (1) 



L{x) = V exp j^'^ Ao(t, x) dt^ 



(2) 



where AJx) is the gauge potential, V denotes path ordering, and r = \x — y\ 



Note that, since center elements of the gauge group commute, the Polyakov loop 
correlator evaluated with center-projected configurations is equal to the Wilson loop 
of identical spatial width and extending along the entire time direction. Therefore, 
center vortices contribute to the static quark potential at finite temperatures in the 
same way as they do at zero temperature (i.e. Stokes' theorem applies in the same 
way). The question to be settled is whether, at finite temperature, vortices still 
provide the entire long-range static quark potential (center dominance) or whether 
other effects become important. This would mean that center vortices lose their role 
as relevant infrared collective degrees of freedom. 

We have tested this empirically employing Monte Carlo measurements. On the one 
hand, we have evaluated Polyakov loop correlators using center-projected links on a 
12^ X Nt lattice with Nt = 5, 6, 7 for different 13 G [2.26 . . . 2.4], i.e. different inverse 
temperatures 1/T = Nta{[3) were achieved by varying (3. The /3-dependence of the 
lattice spacing a was extracted from zero temperature string tension measurements 
(see previous section). While these measurements are thus fraught with sizeable 
uncertainties due to finite size effects in a(/5), the statistical fluctuations still turn 
out to be the dominant source of error. 

The center-projected string tension at as a function of temperature was extracted 
from the Polyakov loop correlators by fitting a linear law to the potential V{r) of eq. 
(p. A Coulomb term is not necessary, since, as will be illustrated by an example 
below, center projection removes the perturbative Coulomb part from V{r) just as 
at zero temperature In addition, however, it should be kept in mind that the 
static quark potential at finite temperatures in general also contains a logarithmic 



dependence on the separation [|T^ . Thus, fitting a purely linear law to the potential 
strictly speaking does not yield the coefficient of the linear term, which by definition 
constitutes the string tension; instead, one obtains an effective "string tension" 
which provides a good parametrization of the long-range static quark potential in 
the limited range of separations accessible to lattice experiments. Since the accuracy 
of our measurements is limited, we cannot separate the linear and logarithmic parts 
of the potential, and we thus quote instead the effective string tension in the sense 
explained above. This quantity in full Yang-Mills theory is known to behave as 
follows. It retains its zero-temperature value to within approximately 10% up to 
the temperature 0.8 T^, where is the deconfinement phase transition temperature, 
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Figure 2: Left: The string tension at calculated from center-projected links as func- 
tion of the temperature T; the four lowest-temperature points on the 12^ x 5 lattice 
are measured at P < 2.3 and thus may already be subject to systematic scaling 
violations. The zero temperature measurement is included to give an indication of 
the error bars. Right: The static quark potential at two sample temperatures, as 
extracted from Polyakov loop correlators with center-projected links, in units of the 
zero-temperature string tension. 



and then quickly drops to zero. For the SU(2) theory, we have only been able to 
find some rather sparse older data to substantiate this |19[; on the other hand, for 



the SU(3) case, which should behave qualitatively in the same manner, new high 
precision measurements exist which corroborate the aforementioned behavior of 
the effective string tension. 

Our results for the center-projected case are displayed in the left-hand plot of Figure 
0. We find that the center-projected string tension at within the error bars repro- 
duces the behavior of the full string tension quoted above; thus, we observe center 
dominance of the long-range part of the static quark potential at finite temperatures 
T < Tc within the accuracy of our measurement and the limited range of separations 
available. 
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Moreover, the center-projected string tension at signals the deconfinement phase 
transition to occur at = 295 ± 10 MeV, where the input scale was a zero tem- 
perature string tension of cr = (440 MeV)^. This value of Tc is in good agreement 
with the high precision measurements of Thus, also the transition to the 

deconfined phase with a vanishing string tension is accurately reproduced by the 
center-projected theory. 

For further illustration, we also calculated Polyakov loop correlators using in par- 
ticular p = 2.3 for a 12^ x A^^ lattice and p = 2.4069 for a 16^ x A^^ lattice while 
varying the number of lattice points Nt in time direction in order to obtain different 
inverse temperatures 1/T = Nta{(3). The j3 values and lattice sizes in both cases 
correspond to a physical size of the spatial cube of Lg = Nsa{(3) ^1.9 fm, while the 
lattice spacing changes by a factor a(/3 = 2.3)/a(/? = 2.4069) = 4/3. This particular 
choice of lattice sizes and /3-values therefore allows to test the scaling behavior of 
the lattice observables with vastly diminished fluctuations due to finite size effects. 
Of course, the overall uncertainty in the mass scale remains. In the case of Polyakov 
loop correlators, this does not turn out to be crucial; the statistical fluctuations 
are the main source of error. However, when measuring vortex densities (see next 
section), this method substantially improves the observed scaling properties, since 
the statistical errors are not as dominant. 

Using the Monte Carlo results for center-projected Polyakov loop correlators on 
these lattices, the static quark potential as a function of the distance at diverse 
temperatures both below and above the deconfinement transition was extracted. 
Two examples, corresponding to Nt = 3,6 for (3 = 2.3 and Nt = 4,8 for (3 = 2.4069, 
are given in the right-hand plot in Figure ^ For T < T^, the potential rises linearly 
even at small distances; as observed previously for the case of zero temperature 0, 
center projection removes the short range Coulomb interaction. 

Vortex polarisation. Given this success of the vortex picture, the question of 
the nature of the deconfinement transition in this picture poses itself. One of the 
simplest hypotheses would appear to be the following: While the center vortices are 
by construction one lattice spacing thick, they represent smooth configurations in 
the original gauge fields before the gauge fixing and center projection procedure is 
applied, with a physical thickness in the continuum limit Possibly vortices 

running perpendicularly to the time direction are simply too thick to fit into the 
space-time manifold of time extension 1/T for T > Tc. This would mean that 
the planar density of points at which vortices pierce the area spanned by the two 
Polyakov loops entering the loop correlator vanishes. Vanishing density of such 
points precludes fluctuations in the number of such points, making an area law decay 
of the Polyakov loop correlator impossible. In order to test this scenario, we have 
measured the planar density pt of points at which vortices pierce an area extending 



8 



a. 



o p_, p=2.3 
- o p,, P=2.4069 
+ P>, P=2.3 
X p,, 13=2.4069 



200 



400 600 

T [MeV] 



800 1000 



1.4 



1.2 



0.8 



0.6 



0.4 



0.2 



o P=2.3, 12 
o p=2.4069, 16' 



200 400 600 800 1000 

T [MeV] 



Figure 3: The planar densities pt and ps of points at which vortices pierce time-hke 
and space-hke planes, respectively. 



in the time and one space direction, for different temperatures. In order to reduce 
fluctuations due to finite size effects (see previous section), these measurements were 
done using (3 = 2.3 for a 12^ x A^^^ lattice and (3 = 2.4069 for a 16^ x A^^ lattice while 
varying the number of lattice points A^^ in time direction in order to obtain different 
inverse temperatures 1/T = Nta{f3). The result for the vortex density is shown in 
Figure ^. 

Evidently, while pt experiences a drop as one increases the temperature past T^, the 
behavior of pt is smooth and at T 2Tc, pt still retains roughly a third of its zero 
temperature value. Thus, the hypothesis advanced above of deconfinement being due 
to a vanishing of pt was too simplistic. There still exists a density of vortices piercing 
the area spanned by two Polyakov loops above Tc, but the random character of the 
distribution of these vortices must disappear in order to realize deconfinement. We 
will come back to this presently; before doing so, note that we have for comparison 
also measured the density ps of vortices piercing an area extending in two spatial 
directions. The ratio pt/ps is also shown in figure ^. For low temperatures, pt and 
Ps coincide, as they must do due to Euclidean 0(4) invariance. At temperatures 
slightly above T^, ps decreases along with pf. At higher temperatures, however, ps 
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begins to increase. This seems consistent with a simple picture of the intersection 
points making up the density ps still being distributed randomly, which leads to a 
linear relation |^ between the planar density ps and the corresponding spatial "string 
tension" as extracted from spatial Wilson loops. This spatial string tension in turn is 
known to increase as ^Jal ~ g'^{T) T for T > 2Tc according to dimensional reduction 
arguments p2[ and their verification in lattice experiments pi| . As mentioned above, 



in the case of the intersection points making up the density pt, this simple random 
picture must by contrast become invalid above Tc. 

Pairing of vortex intersection points. Coming back to the issue of the vortex 
description of the deconfinement transition, it is necessary to inspect more closely 
the properties of vortices piercing the area spanned by the two Polyakov loops en- 
tering a loop correlator. A necessary condition for an area law suppression of the 
loop correlator expectation value is an ever better mutual cancellation, as the area 
spanned by the Polyakov loops is increased, between configurations where the area 
is pierced an even or an odd number of times by vortices, respectively. We have 
therefore measured the probabilities of these two function of tempera- 

ture for an area of spatial width 0.9 fm. For jS = 2.3, this corresponds to 6 lattice 
spacings in the 12^ x A^^ lattice, whereas for (3 = 2.4069, it corresponds to 8 lattice 
spacings in the 16^ x A^^ lattice, i.e. in both cases the distance between the Polyakov 
loops is half the linear extension of the universe. The result of these Monte Carlo 
experiments is shown in Figure^, which displays the fraction p{T) of cases where an 
area specified as above was pierced an even number of times by vortices. This quan- 
tity exhibits a sharp transition at the deconfinement temperature Tc ~ 295 MeV. 
Note that p 7^ 1/2 for large areas precludes an area law decay of the Polyakov loop 
correlator, implying deconfinement. 

The clear signal of the deconfinement transition exhibited by the above probability 
p is not surprising, since it constitutes nothing but a slightly more physical variant 
of the usual Polyakov loop order parameter. The center-projected Polyakov loop 
correlator is given as follows in terms of the probability p that an even number of 
vortices pierces the area spanned by the two Polyakov loops in question: 

L{y)) = p . 1 + (1 - p) . (-1) = 2p-l. (3) 

Thus, p corresponds, up to a rescaling and a shift, to the Polyakov loop correlator 
itself. The only difference to the usual Polyakov loop order parameter, which is 
related to the free energy associated with placing a static color source in the thermal 
medium, is that the static color source has been given a partner (separated from it 
by half the extension of the universe). 

Therefore, in defining the probability p, we have not really introduced a new, dis- 
tinct, order parameter for the deconfinement transition. Rather, the vortex picture 
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Figure 4: Fraction of cases in which an area spanned by two Polyakov loops was 
pierced an even number of times by vortices. 



simply provides an alternative language for describing the Polyakov loop correlator, 
by characterizing the manner in which vortex intersection points occur on the area 
spanned by the two Polyakov loops. The deconfinement transition occurs when the 
intersection points of vortices with the area spanned by the Polyakov loops begin 
to occur predominantly in pairs of finite separation (above Tc), whereas they occur 
randomly below Tc. 

Summary and Outlook. In the present work, we have explored the center vortex 
content of Yang- Mills theory at finite temperatures using numerical lattice methods 
in the direct maximal center gauge and center projection. We have observed vortex 
dominance for the long-renge part of the static quark potential. This establishes the 
relevance of center vortices at finite temperatures and thus provides the necessary 
basis for further investigation of these degrees of freedom. 

The planar densities of intersection points of vortices with a given plane in space- 
time exhibit an anisotropic behavior as temperature is raised from the 0(4)-invariant 
zero temperature limit. Around the deconfinement temperature T^, the density pt 
of intersection points with a plane running in the time and one space direction be- 
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gins to decrease and retains roughly a third of its zero temperature magnitude at 
T = 2Tc. The density pt is relevant for the behavior of the Polyakov loop correla- 
tor, since the area between two Polyakov loops extends in the time and one space 
direction. However, since pt does not vanish for T > Tc, deconfinement cannot be 
simply ascribed to the behavior of pt, i.e. the partial polarization of the vortices in 
time direction. Instead, more detailed correlations between the intersection points 
building up the density pt must be present. 

By contrast, the behavior of the density ps of intersection points with planes running 
in two spatial directions largely parallels the known behavior of the spatial string 
tension ag After a slight dip around Tc, which may be a finite-size artefact. 



Ps begins to rise with temperature along with the spatial string tension as- This is 
consistent with a simple picture of a random distribution of the intersection points 
making up the density ps- In such a picture, the associated (spatial) string tension 
is proportional to the density ||^. 

Returning to the question of the deconfinement transition, we measured in more 
detail the probability that an even or an odd number of intersection points occurs 
on the area spanned by two Polyakov loops at a large distance from one another 
(half the lattice universe). This quantity exhibits the deconfinement transition very 
clearly, which is not surprising, since (up to a trivial shift) it corresponds to the value 
of the Polyakov loop correlator itself, cf. eq. (0). This is a slightly more physical 
realization of the usual Polyakov loop order parameter, namely, instead of putting 
one static color source on the lattice and measuring the associated free energy, we 
give the color source a static partner at the edge of the universe. Thus, we have not 
really defined a new order parameter, but the vortex picture provides an alternative 
language for describing the Polyakov loop correlator: The deconfinement transition 
occurs when the intersection points of vortices with the area spanned by the Polyakov 
loops begin to occur predominantly in pairs of finite separation (above Tc), whereas 
they occur randomly below Tc. We conjecture this behavior to be connected with 
a percolation transition in the vortex surfaces: Since the vortex surfaces are closed, 
intersection points of these surfaces with the entire plane containing two Polyakov 
loops must always occur in pairs. What determines the behavior of the Polyakov 
loop correlator is whether the members of a pair are only separated by a distance 
small compared with the separation of the Polyakov loops, such that most pairs lie 
either entirely inside or outside the area spanned by the Polyakov loops, or whether 
the members of a pair can occur arbitrarily far apart, such that the distribution of 
points is essentially random. In terms of the vortex surfaces, this translates into 
the question whether typical vortex clusters have a finite extension, or whether they 
form networks extending over the entire universe, i.e. percolate. This leads to 
the conjecture that the confining phase is a phase in which the vortices percolate 
throughout space-time, while in the deconfining phase, they cease to do so. We 
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plan to report on a detailed examination of the vortex percolation properties in an 
upcoming publication. 
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